# Generate a single ray in k-space
import numpy as np
from numpy import pi, sin, cos

def inputparsing():
    from optparse import OptionParser

    parser = OptionParser()

    parser.add_option("--nk",action="store",type="int",default="25",help="[%default]")
    parser.add_option("--k",action="store",type="float",default="0.1",help="[%default]")
    op,args=parser.parse_args()

    return op


if __name__=='__main__':

    op = inputparsing()
    k = np.linspace(0.01,op.k,op.nk)


    k  = np.vstack((k.flatten(),k.flatten(),k.flatten()*0+0.05)).T

    np.savetxt('kfile.dat',k)
    
    
